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Context. T Tauri stars exhibit variability on all timescales, whose origin is still debated. 

Aims. We investigate the long term variability of CTTs over up to 20 years, characterize it from a set of statistical parameters 
and discuss its origin. 

Methods. We present a unique, homogeneous database of photometric measurements for Classical T Tauri stars extending up 
to 20 years. The database contains more than 21000 UBVR observations of 72 CTTs. All the data were collected within the 
framework of the ROTOR-program at Mount Maidanak Observatory (Uzbekistan) and together they constitute the longest 
homogeneous, accurate record of TTS variability ever assembled. We characterize the long term photometric variations of 49 
Q ■ CTTs with sufficient data to allow a robust statistical analysis and propose an empirical classification scheme. 
3— j Results. Several patterns of long term photometric variability are identified. The most common pattern, exhibited by a group of 
C/3 , 15 stars which includes T Tau itself, consists of low level variability (AV <0.4mag) with no significant changes occurring from 
season to season over many years. A related subgroup of 22 stars exhibits a similar stable long term variability pattern, though 
with larger amplitudes (up to AV ~1.6 mag). Besides these representative groups, we identify three smaller groups of 3-5 stars 
• i-H , each which have distinctive photometric properties. 

' Conclusions. The long term variability of most CTTs is fairly stable and merely reflects shorter term variability due to cold 
, and hot surface spots. Only a small fraction of CTTs undergo significant brightness changes on the long term (months, years) , 
which probably arise from slowly varying circumstellar extinction. 

Key words, stars: activity of stars - pre-main-sequence stars - variables: general 



1. Introduction 

T Tauri stars ( TTs) wer e first identi fied as a class of ir- 
regular variables by Ijovl lll942lll945h . The properties of 
thes e low mass pre-main seque nce stars were reviewed by, 
e.g., iMenard fc Bertoutl <|l999) who emphasized two sub- 
groups, the so-called "classical" T Tauri stars (CTTs) still 
actively accreting from their circumstellar disks and the 
"weak-line" T Tauri stars (WTTSs) no long er surrounded 
by a circumstellar disk. As discussed by iHerbst et alJ 
l)l994j) . three types of day to week timescale variability 
can be identified in TTs. Type I variability, most often 
seen in WTTSs, is characterized by a low level periodic 
modulation of the stellar flux (a few 0.1 mag) and results 
from the rotation of a cool spotted photosphere. Type II 

* UBVR photometric catalogues are available in electronic 
form at CDS Strasbourg via anonymous ftp to cdsarc.u- 
strasbg.fr. 



variations have larger photometric amplitudes (up to 2 
or 3 mag in extreme cases), are most often irregular but 
sometimes periodic, and are associated with short-lived 
accretion-related hot spots at the stellar surface of CTTs. 
Finally, the more rarely observed Type III variations are 
characterized by luminosity dips lasting from a few days 
up to several months, which presumably result from cir- 
cumstellar dust obscuration, and are associated with the 
so-called UXOrs. 

The ROTOR program (Research Of Traces Of 
Rotation) is dedicated to the study of the photomet- 
ric variability of pre-main-sequence (PMS) and main- 
sequence (MS) objects. The goal of the program is to 
obtain a long term homogeneous database of UBVR 
observations of classical T Tauri stars (CTTs), weak- 
line T Tauri stars (WTTs), post-T Tauri stars (PTTs), 
Herbig AeBe stars (HAeBe), FU Orionis objects, and 
more evolved stars such as RS CVn, BY Dra, and FK 
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Com sta rs. A detailed desc ription of the program was 
given bv IShevchenkol l)l989j) . More than 100000 obser- 
vations of about 370 variable stars have been obtained 
since 1983. Our list contains 72 CTTs, 45 WTTs, 54 
PTTs, 88 HAeBe, and 11 FK Com stars. Observations 
were secured over more than 20 years by members of 
the Tashkent Astronomical Institute, including S.Yu. 
Melnikov (24 % of all observations), K.N. Grankin (23 %), 
M. Ibrahimov (14 %), S.D. Yakubov (10 %), O.V. Ezhkova 
(10 %), and V. Kondratiev (7 %) . Preliminary results were 
published bv IShevchenko et all lll993blal). iGrankin. et alJ 
1995HShevchenko et al 111 998UOra,nkinl ||1 9971 ll 99Slll 999h . 
Herbst fc Shevche nko (199.9J). 

In this paper, we analyze the unique and extensive 
Mt. Maidanak database to investigate CTTs variability 
over a timescale of several years and, in some cases, sev- 
eral decades. In Section 2, we describe the stellar sample 
and the observations carried out for the last 20 years at 
Mt. Maidanak. In Section 3, we define the statistical pa- 
rameters used to characterize the long term variability of 
CTTs. In Section 4, we provide the results of the statisti- 
cal analysis of CTTs light curves and offer in Section 5 an 
empirical classification of their long term photometric pat- 
terns. In Section 6, we attempt to relate the various types 
of photometric patterns seen to the underlying causes of 
variability. 

2. Star sample and observations 

All UBVR data were obtained at Mount Maidanak 
Observatory (longitude: E4 h 27 m 35 s ; latitude: +38°41'; al- 
titude: 2540 m) in Uzbekistan. Results of astroclimate 
st udies at Maidanak observat ory have been summarized 
bv lEhgamberdiev et alJ 1 2000l) . who show that it is among 
the best observatories in the world, in many respects. We 
targ eted 72 CTT s for lo ng-term monitoring from the list 
of lHerbig fe Bell (|l988j) . More than 21000 UBVR mag- 
nitudes were collected for these objects from 1983 up to 
2003, although the number of U magnitudes is relatively 
small compared to the other colors due to the faintness of 
the stars at those wavelengths. All the observations were 
obtained at two telescopes (0.6 and 0.48 m reflectors) us- 
ing a single-channel pulse-counting photometer with pho- 
tomultiplier tubes. As a rule, each program star was mea- 
sured once per night at minimal airmass. We secured up 
to 120 (U)BVR measurements for each sample target dur- 
ing each observational season lasting several months. The 
number of measurements depends on the object's visibility 
from Mt. Maidanak. The largest number of measurements 
was thus obtained for CTTs in Cygnus, Cassiopea, and 
Taurus, and fewer in Ophiuchus and Orion. Magnitudes 
are provided in the Johnson UBVR system. The rms error 
of a single measurement in the instrumental system for a 
star brighter than 12™ in V is about m .01 in BVR and 
m .05 in U. 

Observations were carried o ut either differ entially us- 
ing a nearby reference star (see IStrauzislll977l) or dire ctly 
by estimating the nig htly extinction l|Nikonovl It976|) . In 



the latter case, several reference stars were observed ev- 
ery night to derive the exti nction coefficients in each fil- 
ter. Selected standard stars l|Landoltlll983Lll992() were ob- 
served and used to calibrate instrumental magnitudes on 
the Cousins system. We then transformed the magnitudes 
to the Johnson U BVR system using the relationship from 
lLandoltl lll983l) : (V-R) c = -0.0320 + 0.71652x (V-R) j. 
The formal accuracy of this reduction step is m .01. All 
observed local times have been converted to heliocentric 
julian days (JDH). A detailed description of the equip- 
m ent, observat i on tec hniques, and data processing is given 
in IShevchenkol l|l989lh 

The results of the Maidanak program of homoge- 
neous long-term photometry of CTTs are summarized 
in Table 1. Column s are: star's name, star's number in 
iHerbig fe Belli l|l988j) 's catalogue, spectral type, time span 
of observations in JD, number of seasons observed, photo- 
metric range in the V band, number of observations in the 
V band, and average values of B — V and V — R colors. 
The complete Maidanak database for CTTs is available in 
electronic form at CDS, Strasbourg, via anonymous ftp to 
cdsarc.u-strasbg.fr. 

3. Statistical parameters 

We define here the statistical parameters which allow us 
to characterize the long term photometric variability of 
CTTs. To increase the reliability of the statistical analy- 
sis we consider only observing seasons having at least 12 
measurements and further exclude those stars which have 
less than 5 seasons of observations. As a result we retain 
49 CTTs for the analysis. We also performed a more strin- 
gent analysis by retaining only 39 stars having more than 
20 measurements per season and more than 5 seasons of 
observations. The results of the 2 analyzes are not signifi- 
cantly different from one another. We thus conclude that 
our selection criteria for the first sample of 49 stars provide 
robust statistical parameters (mean, standard deviations, 
minimum and maximum levels, etc.). 

For each object, we compute the mean light level 
(V m ^') in season (i) and the corresponding photomet- 
ric range (AVW) for that season. Then we average the 
(V m ®) over all seasons to yield (V m ) and also compute 
its standard deviation (ov m ) based on the scatter of the 
seasonal means. Similarly, we compute the average V band 

Y AV (i) 

photometric amplitude over all seasons (AV — N , 

where N s is the number of observational seasons) and its 
standard deviation [aw)- We also use the ratio ^pr- 

In addition to statistical parameters, we investigate 
the color behavior of CTTs. Previous studies of the short 
term variability of CTTs have shown that most stars tend 
to redden in color indices (such as V — R and B — V) when 
fading (e.g. iHerbst et al.| [l9941 This beha vior can result 
from surface spots, either cold or hot (cf. iBouvier et alJ 
1199.1) . or from variable extinction by circumstellar dust. 
We characterize the color behavior of CTTs by comput- 
ing the color slopes A ^y V ^ and A ^y R ^ f rom least-square 
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Fig. 1. The long term V-banel light curve of the CTTs UY 
Aur. Top panel: photometric measurements; lower panel: 
statistical analysis (see text for details). 

fits, and the correlation coefficients pB-v an d Pv-R- We 
used the full set of data pertaining to each object in or- 
der to compute these parameters since there are usually 
too few measurements per season to obtain accurate sea- 
sonal color slopes. For most stars, a strong correlation is 
found between color and brightness variations with little 
evidence for significant variations from season to season 
(cf., e.g., Fig. 8). 

The resulting values of statistical and color parame- 
ters for 49 CTTs are listed in Table 2. The columns in 
Table 2 provide the star's name, the number of seasons 
N* with at least 12 measurements, the mean brightness 
level (V m ) and its seasonal standard deviation (ov m ), the 
average photometric range (AV, averaged over all seasons) 
and its seasonal standard deviation (<tav), the fractional 
variation of photometric amplitude as measured by the 
ratio and the color slopes A ^y^ and A ^v^ to- 

gether with their correlation coefficients ps-v an d Pv-R- 
Two additional parameters discussed below (cf. Sect. 4) 
are also listed, which characterize the relative seasonal 
variability of the maximum and minimum light levels 
(C2 = oy min l<Jy max ) and the preferred brightness state of 
the object on the long term (CI =< v ™' d ~v™™ > where 
V me d is the median light level of the object in a given sea- 
son and the bracket indicates an average over all seasons). 
For instance, a low CI value indicates an object which 
spends more time close to its maximum brightness level. 

In order to illustrate the meaning of these parame- 
ters, an example of the statistical analysis is shown in 
Figure 1 for UY Aur. The long-term V-band light curve 
obtained at Maidanak Observatory is shown in the upper 
panel. In the lower panel, some of the statistical parame- 
ters are illustrated. The mean light level for each season 
of observations (V m ^) is marked by a filled circle and the 
minimum and maximum brightness by crosses. Statistical 
parameters are not computed for two seasons (2001 and 
2004) which contain less than 12 measurements. The solid 



line joins values of average brightness and demonstrates 
significant changes on years timescale. The photometric 
amplitude is also seen to vary on the long term. 

Figures 2a and 2b display the light curves of the other 
48 CTTs on which we applied the same statistical anal- 
ysis. It is immediately obvious that different photometric 
patterns occur on the long term. We discuss these patterns 
with the help of the above defined statistical parameters 
in the next section. 

4. Results 

In this section, we present the results of the statistical 
analysis of the long term light curves. The histogram of 
the average V-band amplitude of variability for the CTTs 
in our sample is shown in Figure 3. The distribution is 
rather large, with AV ranging from m .l up to 2 m . 

Nearly half of the CTTs (21) exhibit small ranges of 
variability on the long term, from m .l to m .4. A roughly 
equal number (24) display larger ranges, from 0"\6 to l m .5 
in the V-band. Only a few CTTs (4) exhibit average V- 
band ranges larger than l m .5 and up to 2 m over the years. 
These are BM And, CI Tau, AS 205, and DK Tau. 

A histogram showing the frequency distribution of the 
measured variations of mean light level from season to 
season of the 49 CTTs in our statistical sample is displayed 
in Figure 4. Specifically, this variation is measured by the 
standard deviation of the seasonal mean light levels (cry m ) . 
The histogram is strongly peaked, with the large majority 
of CTTs demonstrating little variation in their mean light 
level over the years (<Jv m — m .l — m .2). However, a 
few CTTs (6) are characterized by significant changes in 
average brightness, with oy m ranging from m .25 to m .65 
(cf. Fig. 6). These are CO Ori, DG Tau, HN Tau, UY Aur, 
RY Tau, and HI Ori. Figures 1 and 2 show that the mean 
light level can vary smoothly on a timescale of several 
years (e.g. UY Aur) or, on the contrary, can change quite 
abruptly from one season to the next (e.g. DG Tau). 

Similarily, the variations of the V-band range from sea- 
son to season are small for most CTTs, amounting to only 
a few tenths of a magnitude at most, as shown by the his- 
togram of (Jav m Figure 5. We caution that in seasons 
where the number of data is not large or in the case of stars 
where the variations are rapid and short-lived, we may be 
underestimating the true photometric range. Nonetheless, 
it appears that this histogram is not as peaked towards 
small values as is the histogram of oy m (Fig. 4) but dis- 
plays a more continuous distribution. We find that the rel- 
ative changes in amplitude (==) are usually constrained 
to lie between 0.1 and 0.5, except for 4 CTTs which do 
exhibit larger seasonal variations up to 0.9 (EZ Ori, SU 
Aur, V2062 Oph, and GW Ori, see Fig.6). In other words, 
CTTs with the largest photometric amplitudes are also 
those whose amplitude varies the most on a timescale of 
years. 

That most CTTs exhibit relatively similar variations 
in mean brightness level and fractional amplitudes is best 
seen in Figure 6 where the 2 quantities are plotted against 
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Fig. 2a Long term V-band light curves of 48 CTTs observed at Mt. Maidanak. 
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Fig. 2b Statistical representation of the 48 CTTs light curves shown in Fig. 2a 
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Fig. 3. Average photometric amplitude in the l/-band for the Fig. 4. Standard deviation of the seasonal mean brightness 
49 CTTs of our sample (M 7 ). level {ov m ). 
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Fig. 5. Standard deviation of the seasonal photometric Fig. 6. The relative variation of amplitude vs. the variation 



amplitude (jav). 



of the mean brightness level. Four objects with the largest 
relative variations of amplitudes are shown by empty circles. 
Five stars with the largest changes in mean light level are 
denoted by asterisks. 



each other. The vast majority of CTTs in our sample lie ing more time i n the high luminosity state than in the low 

in a small region of the diagram delimited by ay m < m .3 one. IParenaeol l|l954f) discussed a similar classification of 

and Z= < 0.5. The few deviant CTTs, already men- irregular variables based on a histogram of their measured 

tionned above, are easily identified in this plot. brightnesses. 



We define 2 additional parameters in order to char- 
acterize the preferred brightness state of each object and 
the long term variations of its extreme (mimimum and 
maximum) light levels. The first parameter is defined as 
CI =< v ™ed-v min > where V me d is the median bright- 
ness of the object in a given season, V m i n its maximum 
brightness level and AV^ its photometric amplitude dur- 
ing this season. CI is computed for each season and then 
averaged over all seasons. The CI parameter thus charac- 
terizes the preferred light level of the object on the long 
term. Thus defined, low CI values indicate objects spend- 



The second parameter is defined as C2 = oy min / ov max 
where V m in is the maximum brightness state and V max the 
minimum one for each season. This parameter thus charac- 
terizes the long term variability (years) of the bright state 
relative to that of the faint state. The C2 parameter can be 
useful to clarify the origin of the variability. For instance, 
in CTTs where the variability primarily results from hot 
spots, the minimum light level (~ photospheric level) is 
expected to remain nearly constant over the years, while 
the maximum light level (star + hot spot) is expected to 
vary on the long term as the spot properties (size and/or 
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Fig. 7. The star's preferred brightness state (as measured 
by Cl=< V " le '^ /m '" >) is plotted against the relative 
variations of its extreme brightness states (as measured by 
C2=ov mi „/ov„ lax ). See text. 
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Fig. 8. Illustration of the linear relation found for most objects 
between the V — R color and V magnitude. Color slopes and 
correlation factors are listed in Table 2. 
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Fig. 9. Illustration of the non linear V — R color variations on 
V magnitude found for a few objects. 




Fig. 10. Least-square color slopes. Top panel: A ^ V V ^ ; lower 
panel: A( Y~ fl) . 



temperature) change. Hence, CTTs with hot spots are ex- 
pected to have large C2 values. Conversely, CTTs whose 
variability is due to cold spots or circumstellar extinction 
will exhibit a nearly stable maximum light level (~ pho- 
tospheric flux) while the minimum light level is expected 
to vary as the cold spot or the occulting blob size evolve 
over the years. Hence, CTTs whose variability is due to 
cold spots and/or circumstellar extinction are expected to 
have low C2 values. 

CI and C2 were computed for 35 CTTs observed over 
at least 7 seasons, in order to derive a reliable estimate of 
<J v min /^v max ■ Figure 7 shows CI against C2 and reveals 



some correlated behavior of the 2 parameters. The group 
of objects in the lower left part of the diagram spend most 
of their time close to the bright state (low CI). They also 
exhibit a fairly constant maximum light level, while large 
variations of minimum brightness occur on the long term 
(see, e.g. AA Tau and GW Ori in Fig. 2). This behavior is 
suggestive of non uniform circumstellar extinction occur- 
ring on the short term (weeks) with an amplitude which 
varies on the long term (years). In support of this hypoth- 
esis, we note that 8 of 10 UX Ori objects monitored at 
Mt Maidanak lie in the same region of the diagram (not 
shown). Also, the low frequency of such objects (5/35) 
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agrees well with the probability of observing oc cultation 
events from the circumstellar disk as computed bv lBertoutl 
(200(3). Alternatively, slowly evolving cold surface spots 
may produce this pattern for highly inclined objects. For 
this group of objects, the maximum light level corresponds 
to the photospheric luminosity. 

At the other extreme, objects in the upper right part 
of the diagram spend most of their time close to the faint 
state (high CI) and their maximum light level varies sig- 
nificantly on the long term (see, e.g. DF Tau and DN 
Tau in Fig. 2). The variability of these objects is probably 
dominated by hot surface spots whose properties vary over 
the years. The photospheric luminosity of these objects 
would thus be best measured close to minimum brightness. 
Finally, most CTTs lie in the intermediate region of the 
diagram. They do not have any preferred brightness state 
(Cl~0.5) and their extreme brightness levels exhibit the 
same amount of long term variability (C2~l), and some- 
times (but not always) vary in a correlated fashion. The 
photometric variability of these objects may result from a 
mixture of hot and cold surface spots whose surface distri- 
bution and/or properties vary on the long term (see, e.g. 
DR Tau, T Tau, BM And, DI Cep in Fig. 2). We show in 
an accompanying paper that WTTS, whose photometric 
variability is due to slowly evolving cold surface spots, lie 
in this intermediate region of the diagram as well. 

Finally, we investigate the color behavior of the 49 
CTTs in our sample. For most objects, the colors (ei- 
ther (V — R) or (B — V)) scale linearly with brightness. 
Hence, the color slope can be computed from a linear least 
square fit. For most CTTs (39), the correlation coefficient 
between V — R and V variations is higher than 0.47 (cf. 
Fig. 8). Only 5 stars do not show any reliable dependence 
of V — R on brightness (correlation factor of lower than 
0.47) when averaged over all seasons. These stars are: GG 
Tau, LkHa 191, LkHa 321, LX Ori, and V360 Ori. In ad- 
dition, 5 CTTs demonstrate a more complex color behav- 
ior, some exhibiting a blue turnaround towards minimum 
brightness (cf. Fig. 9). The unique color pattern of these 
stars (V521 Cyg, RY Tau, EH Cep, AA Tau, and CO Ori) 
is discussed below. 

The histogram of A ^ y ^ color slopes for the 39 CTTs 
with linear color variations is shown in Figure 10. The 
histogram peaks at values of order of 0.2-0.3, close to the 
value expected for variable circumstellar extinction (0.25), 
assuming interstellar extinction properties. This could be 
taken as evidence that CTTs light variability is caused by 
a changing screen of circu mstellar dust i ntercepting the 
line of sight to the star fe.g. lBertoudl200o() . However, sim- 
ilar color slopes can be reproduced as well by other mech- 
anisms. We show in an accompanying paper that WTTS 
exhibit the same average A ^y R ^ color slope as CTTs. 
In WTTS, brightness and color variations result from the 
rotational modulation of the stellar flux by cold surface 
spots. Hence, color slopes usually do not uniquely distin- 
guish between several possible sources of variability. 



A similar analysis was carried out for the A ^ color 
slope, whose histogram is also shown in Figure 10 and 
looks quite similar to the redder color, with a maximum at 
0.2. The ISM extinction law for A ^v^ predicts a value 
of 0.32 if color changes are due to variable circumstel- 
lar extinction. However, besides extinction, other sources 
may affect the bluer colors of CTTs, for example the well 
known veiling phenomenon which is most important in the 
blue. Here again, color slopes are not sufficient to distin- 
guish between the various possible sources of photometric 
variations. 



5. Discussion 

The short term (weeks) photometric variability of CTTs 
is usually interpreted as resulting from 3 possible causes, 
sometimes present simultaneously: (i) cold surface spots, 
presumably resulting from strong photospheric magnetic 
fields; (ii) hot spots at the stellar surface, related to disk 
accretion onto the star; and (iii) variable obscuration of 
the star by circumstellar matter. 

The results obtained from the statistical analysis of 
the long term light curves of 49 CTTs are graphically 
summarized in the 3D plot shown in Fig. 11. From these 
results, we attempt to identify sub-groups of stars which 
appear to share the same type of long term photometric 
behavior and try to relate it to specific physical processes. 
A major difficulty is that different processes may be at 
work on different timescales (hours, days, weeks, years), 
and several sources of variability are likely to coexist. In 
addition, other parameters, such as the system inclina- 
tion, may yield a different photometric pattern for a given 
source of variability. We thus merely try to identify the 
dominant source of variability for each subgroup, among 
the 49 stars with sufficient data. 

The largest group of CTTs we identify in Fig. 11 con- 
sists of 37 stars which exhibit homogeneous long term 
photometric variability. Of these, 15 are located close to 
the origin of the plot exhibiting a low level of variability. 
These 15 stars are shown by filled squares in Fig. 11 and 
the other 22 stars as crosses. The spectral types of these 
stars span a range from G5 up to M0. This group includes, 
e.g., T Tau and DI Cep. CTTs in this group are character- 
ized by relatively stable patterns of variability on the long 
term. Many are known to be actively accreting from their 
disk and the analysis of their short term (weeks) photo- 
metric variability suggests it is dominated by hot surface 
spots, though cold s pots as well as circumstellar extinction 
can also be present llBouvier fc BertoulJIl 989t I Vrha et al] 
Il993t iBouvier et alJll99.4 lBertoutJl2000ji . Hence, the long 
term variability of these stars merely reflects variability 
on shorter timescales, with no significant variations from 
season to season. 

The constancy of the light curves of these objects over 
the years suggests that even though the properties of sur- 
face hot spots appear to vary on the short term, possibly 
due to the changing accretion rate onto the star, the accre- 
tion process remains active and fairly stable on a timescale 
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Fig. 11. A 3D plot of the statistical properties of the long term light variability of 49 CTTs. The compact group of stars located 
close to the origin (filled squares) exhibit the lowest degree of variability. Five stars demonstrate large changes in their mean 
brightness level over the years (asterisks). Another small group exhibit significant seasonal variations in their photometric 
amplitude (circles). Five CTTs exhibiting unusual color behavior with a blue turnaround at minimum brightness are shown as 
triangles. 



of several years. The large range of photometric amplitude 
displayed by CTTs in this group, from a m .l up to l m 6 
in the V-band, may reflect either different accretion rates 
onto the star or geometric projection effects resulting from 
the combination of the spots latitude and the system's in- 
clination. We do not find any correlation between photo- 
metric amplitudes and the time-averaged accretion rate in 
this sample using the H Q flux as a proxy, which favors the 
latter explanation. 

Besides this main group, a much smaller subgroup of 
deviant CTTs seen in Fig. 11 exhibits large changes in 
their mean brightness level over the years (DG Tau, HN 
Tau, UY Aur, CO Ori, and RY Tau). These stars are 
shown by asterisks in Fig. 11. Some of these, like UY Aur, 
display a smooth change of mean brightness over the years 
(Fig. 1). Others however, exhibit abrupt mean brightness 
variations from one season to the next (e.g. from ll m .71 
to 13 m .31 between 1997 and 1998 for DG Tau). The ori- 
gin of these variations cannot be unambiguously assigned 
to a single mechanism. Large scale variations of the ac- 



cretion rate onto the star could produce such a behavior. 
Alternatively, a variation of the circumstellar extinction 
on the line of sight on a timescale of years cannot be ex- 
cluded. Such a timescale would then point to an obscuring 
screen located in the outer disk regions. 

Another small group exhibits significant seasonal vari- 
ations in their photometric amplitude (EZ Ori, SU Aur, 
V2062 Oph, and GW Ori). We have shown these stars as 
circles in Fig. 11. The long term light curve of these objects 
is characterized by a nearly constant maximum brightness 
level with a usually small amplitude of variability, but in- 
terrupted at times by deep fading episodes. The amplitude 
changes from season to season can amount to 52-90% of 
the average amplitude in these objects, while it amounts 
to only 10-47% in the rest of the CTTs sample. Some of 
these objects might be eclipsing binaries (e.g. GW Ori). 
Others may undergo partial occultation by circumstellar 
disk material on a timescale ranging from weeks to years 
(e.g. SU Aur, cf. IPeWarf et~aT]l2008L lBertoutJl2000l) . 
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Finally, a last group of a few CTTs exhibit unusual 
color behavior with a blue turnaround at minimum bright- 
ness (V521 Cyg, RY Tau, EH Cep, AA Tau, and CO Ori). 
These stars are shown as triangles in Fig. 11. This likely 
points to variable circumstellar extinction : as the stellar 
photosphere is partly occulted by circumstellar material, 
the ratio of scattered light to direct light increases, and 
the system turns bluer. These stars tend to be among the 
earlier type stars in our sample, with spectral types from 
F8 to K2, except f or AA Tau, quite a typica l CTT s with 
a spectral type K7 l|Bouvier et alJll999l I2OO31 12006^ . Such 
occultation effects are expected to be observed i n only a 
small fraction of CTTs seen at high inclination I Bertoutl 
l2000h . Two stars of this group (RY Tau and CO Ori) also 
belong to the small group of CTTs with large changes in 
their mean brightness level over the years. They are shown 
by two superimposed symbols (asterisks and triangles) in 
Fig. 11. 

Our 20-year long photometric database thus provides 
the first opportunity to investigate the long term variabil- 
ity of a significant sample of CTTs. Yet, this baseline re- 
mains short compared to the accretion (viscous) timescale 
in an accretion disk, not to mention the Kelvin-Helmoltz 
PMS evolution timescale. Hence, while the longest pho- 
tometric database available today on PMS stars, it may 
still miss some sources of variability acting on even longer 
timescales. As an illustr ation, the 40-yr long l i ght c urve 
of T Tau investigated bv lMernikov fc Grankml pOOfil) ex- 
hibits a gradual rise of the mean light level from V~10.5 
to V~9.9 over a timescale of about 20 years prior to 1980. 

6. Conclusion 

The photometric variations of a large sample of CTTs 
have been monitored at Mt. Maidanak Observatory. More 
than 21 000 UBVR photometric measurements have been 
collected in a uniform manner for 72 CTTs over as many 
as twenty years. Based on this unique database, we in- 
vestigated here the long term photometric properties of a 
subsample of 49 CTTs that were observed at least 12 times 
in each of 5 or more seasons. By defining a set of statistical 
parameters, we characterized the long term photometric 
patterns of this sample of CTTs. 

Several types of long term photometric patterns were 
identified. The most common type, shared by about 75% 
of the CTTs in our sample, is characterized by low to 
moderate amplitudes of variability and relatively stable 
light variations on a timescale of several years. We ar- 
gue that the source of variability in these objects is likely 
dominated by a mixture of hot and cold surface spots. 
In this subgr o up, 15 stars were previously investigated by 
iHerbst et all l|l994|) and classified as Type II variables. 
Even though the properties of the spots may change on the 
short term (weeks), the overall stability of the light vari- 
ations over the years suggest that the underlying physical 
processes, i.e., accretion onto the star for hot spots and 
strong magnetic fields for cold spots, remain instrumental 
and undergo only marginal changes on the longer term. 



In addition to this common pattern, a few more pecu- 
liar types of long term photometric variations are found. 
They are characterized by significant changes in the mean 
brightness level and/or photometric amplitudes over the 
years. We suggest that, in most cases, this type of un- 
usual variability is dominated by variable circumstellar 
extinction, presumably due to the inner (short term) or 
outer (long term) disk regions partially occulting the cen- 
tral star. Occultation effects are likely to happen only in 
objects seen at high inclination, which e xplains why the se 
patterns are more rarely observed (cf. iBertoutl (2000)), 
being exhibited by only 12 CTTs over 49 in our sample. 
Most of the stars in this s ubgroup correspond to Type 
III variability as defined bv lHerbst et alJ (|l994). As pre- 
viously noted by these authors, we find that most (8/12) 
of the CTTs exhibiting this peculiar type of variability 
have a spectral type K2 or earlier, while most stars in our 
sample have later spectral types. Hence, this type of vari- 
ability seems to preferentially (but not exclusively) affect 
early-type CTTs and may be related to the UXor variabil- 
ity pattern seen in a minority of the more massive Herbig 
Ae-Be stars. 

Finally, we caution that photometric variability in 
CTTs, be it short or long term, often results from sev- 
eral processes acting simultaneously and which may be 
difficult to identify from photometry alone. Spectroscopic 
and polarimetric measurements obtained simultaneously 
with photometric ones, which have been gathered so far 
for very few CTTs, usually offer a clearer view of the phys- 
ical processes responsible for the observed variability. 

The results reported here on the long term photomet- 
ric behavior of CTTs will be contrasted with a similar 
analysis of the long term photometric variations of WTTs 
from the same database in an accompanying paper. 
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Table 1. Results of Maidanak long-term photometry of CTTs 



Name 


HBC 


SpT 




N s 


V range 


N obs 


B — V 


V-R 


AA Ori 


130 


K4 


47055-47777 


1 


12.79-13.16 


11 


1 


.28 


1 


.23 


AA Tau 


63 


K7 


46690-53028 


15 


12.30-14.14 


568 


1 


.34 


1 


.40 


AS 205 


254 


K5 


46955-49588 


8 


10.87-13.77 


362 


1 


.24 


1 


.55 


AV Ori 


159 


K4 


47040-47173 


1 


12.84-14.42 


14 


1 


.23 


1 


.23 


BD-22 4059 


608 


G5 


47670-49514 


5 


10.10-10.28 


165 





.84 





.75 


BM And 


318 


K5 


47008-52908 


17 


11.85-13.98 


999 


1 


.13 


1 


.07 


BP Tau 


32 


K7 


46679-52938 


11 


11.67-12.99 


495 


1 


.06 


1 


.26 


CI Tau 


61 


K7 


47018-52935 


7 


12.28-13.79 


319 


1 


.37 


1 


.55 


CO Ori 


84 


F8 


47031-52927 


16 


9.81-12.73 


541 


1 


.06 





.99 


CX Tau 


27 


M2.5 


51782-51841 


1 


13.62-13.80 


15 


1 


.51 


1 


.65 


CY Tau 


28 


Ml 


52579-53028 


1 


12.94-13.59 


19 


1 


.07 


1 


.53 


DD Tau 


30 


Ml 


48123-49275 


4 


13.45-15.15 


114 


1 


.07 


1 


.76 


DE Tau 


33 


M2 


46681-47382 


1 


12.82-13.95 


23 


1 


.43 


1 


.64 


DF Tau 


36 


M0 


45920-53028 


19 


10.64-12.74 


929 


1 


.10 


1 


.45 


DG Tau 


37 


M 


45928-52938 


18 


11.23-13.58 


710 


1 


.06 


1 


.35 


DH Tau 


38 


Mle 


52932-52938 


1 


13.20-13.73 


5 





.89 


1 


.36 


DI Cep 


315 


G8 


46966-52907 


17 


11.14-11.67 


1231 





.89 





.88 


DI Tau 


39 


M0 


46710-51841 


11 


12.69-12.97 


435 


1 


.59 


1 


.50 


DK Tau 


45 


K7 


47022-52937 


11 


11.54-14.46 


429 


1 


.27 


1 


.43 


DL Tau 


58 


K7 


47017-52935 


8 


12.39-13.89 


370 


1 


.09 


1 


.42 


DM Tau 


62 


Ml 


47028-52935 


6 


13.43-14.32 


162 


1 


.07 


1 


.51 


DN Tau 


65 


M0 


46681-52914 


10 


12.08-12.54 


369 


1 


.33 


1 


.35 


DR Tau 


74 


K5 


46686-52924 


16 


10.76-12.78 


642 





.82 


1 


.05 


DS Tau 


75 


K5 


46680-48230 


4 


11.58-12.68 


199 





.91 


1 


.09 


EH Cep 


307 


K2 


46970-52907 


17 


11.52-13.44 


1296 


1 


.19 


1 


.03 


EZ Ori 


114 


GO 


47048-53022 


7 


11.48-12.25 


222 





.84 





.77 


GG Tau 


54 


K7 


46686-49287 


7 


11.91-12.48 


345 


1 


.42 


1 


.44 


GI Tau 


56 


K6 


46681-52935 


8 


12.30-14.85 


344 


1 


.46 


1 


.57 


GK Tau 


57 


K7 


46682-52935 


8 


11.94-13.58 


376 


1 


.43 


1 


.44 


GM Aur 


77 


K3 


46681-50046 


9 


11.74-12.35 


349 


1 


.21 


1 


.22 


GW Ori 


85 


G5 


47031-52927 


15 


9.74-10.53 


531 





.99 





.94 


GX Ori 


89 


Kl 


47031-47169 


1 


13.11-13.47 


20 





.95 


1 


.02 


Haro 1-1 


256 


K5.7 


46964-49190 


6 


12.80-13.82 


166 


1 


.27 


1 


.45 


Haro 1-14 


267 


M0 


46964-50660 


6 


13.75-15.03 


154 


1 


.52 


1 


.68 


Haro 1-4 


257 


K6.7 


46964-48837 


5 


13.29-13.81 


131 


1 


.75 


1 


.85 


Haro 1-8 


261 


K5 


46965-49190 


4 


13.74-14.22 


133 


1 


.69 


1 


.70 


Harol-14c 


644 


K3 


46607-46649 


1 


12.13-13.13 


27 


1 


.67 






HI Ori 


93 


K4 


47030-52176 


7 


12.63-13.86 


214 


1 


.28 


1 


.20 


HL Tau 


49 


K7-M2? 


48147-48954 


2 


14.22-14.63 


48 


1 


.42 


1 


.54 


HN Tau 


60 


K5 


47028-49277 


6 


12.63-14.94 


163 





.97 


1 


.30 


HP Tau 


66 


K3 


46702-47421 


1 


13.35-14.20 


35 


1 


.78 


1 


.69 


LkCa 8 


385 


M0 


48953-52935 


1 


12.83-13.44 


26 


1 


.31 


1 


.39 


LkCa 15 


419 


K5 


48859-50760 


5 


11.85-13.34 


123 


1 


.26 


1 


.16 


LkHa 172 


298 


K4 


48087-49962 


2 


14.58-14.98 


64 


1 


.44 


1 


.36 


LkHa 191 


301 


K0 


46618-49956 


5 


12.87-13.09 


326 


1 


.10 


1 


.04 


LkHa 200 


6 


Kl 


46258-46802 


1 


13.51-13.68 


29 


1 


.23 


1 


.10 


LkHa 209 


194 


G8 


46722-50044 


8 


11.64-12.06 


193 





.83 





.80 


LkHa 228 


295 


K1.2 


46261-50017 


2 


12.87-13.06 


48 





.84 





.78 


LkHa 274 


203 


K4 


46722-50024 


3 


13.16-13.49 


92 


1 


.26 


1 


.20 


LkHa 321 


303 


G 


46985-50020 


7 


12.25-12.44 


547 


1 


.23 


1 


.17 


LX Ori 


133 


K3 


47039-51834 


6 


11.98-12.34 


191 


1 


.10 


1 


.00 


Oph 6 


653 


K 


48419-49245 


1 


13.18-13.45 


48 


1 


.69 


1 


.75 


RW AurAB 


80.81 


K1-K3 


46680-52924 


16 


9.32-11.75 


534 





.72 





.87 


RY Tau 


34 


Kl 


45968-53022 


19 


9.55-11.56 


885 


1 


.03 


1 


.07 
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Table 1. continued. 



Name 


HBC 


SpT 




N s 


V range 


N obs 


B — V 


V — R 


S CrA 


286 


K6 


47673-48838 


3 


10.40-12.32 


122 


1.00 


1.18 


SU Aur 


79 


G2 


45676-52926 


19 


8.92-10.02 


793 


0.92 


0.83 


SW Ori 


115 


K4 


47041-50022 


8 


12.42-13.69 


153 


1.04 


0.97 


T Tau 


35 


K0 


46665-52937 


17 


9.75-10.18 


982 


1.17 


1.13 


UX TauAB 


42.43 


M1-K2 


47378-49287 


6 


10.64-12.77 


318 


1.17 


1.06 


UY Aur 


76 


K7 


46681-52914 


14 


11.28-13.68 


427 


1.27 


1.43 


UZ Tau 


52 


Ml. 3 


46383-47101 


1 


12.25-13.11 


15 


1.19 


1.66 


V1082 Cyg 


728 




48059-50022 


6 


12.62-13.67 


461 


0.76 


1.00 


V1121 (Jpn 


270 


K5 


46955-52876 


14 


1 A O A 11 OO 

10.80-11.83 


a ai 

601 


1.28 


1 on 

1.36 


V2058 Oph 


259 


K6 


46607-49194 


7 


12.60-13.09 


209 


1.49 


1.60 


V2062 Oph 


268 


K2 


46607-49194 


7 


12.20-13.05 


223 


1.34 


1.39 


V360 Ori 


144 


K6 


47060-50022 


5 


12.40-12.71 


159 


1.39 


1.30 


V521 Cyg 


299 


G8 


46609-52899 


17 


13.37-14.70 


1140 


1.27 


1.20 


V625 Ori 


183 


K6 


47061-50010 


5 


12.61-13.97 


114 


1.08 


1.15 


V649 Ori 


86 


G8 


47031-52176 


7 


11.87-12.31 


275 


1.12 


1.09 


V853 Oph 


266 


Ml. 5 


46607-51014 


7 


12.61-13.87 


238 


1.10 


1.63 


XZ Tau 


50 


M3 


48142-48954 


3 


12.96-15.17 


74 


1.40 


1.88 


YY Ori 


119 


K5.M0 


47041-47170 


1 


13.20-14.08 


16 


0.77 


1.01 
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Table 2. Statistical properties of CTTs light curves (see text). 



Star Name 


iv; 




V m 


°v m 


AV 


&AV 


"AV 
AV 


A(B-V) 
AV 


pB-V 


A(V-B.) 
AV 


Pv-b 


cl 


c2 


AA Tau 


15 


12 


917 





143 


1 


133 





322 





285 





017 





04 





093 





54 





423 





304 


AS 205 


8 


12 


259 





224 


2 


002 





459 





229 





129 





51 





239 





89 





997 





452 


BD-224059 


5 


10 


169 





016 





116 





035 





300 





073 





08 





483 





53 





525 





427 


BM And 


17 


12 


809 





189 


1 


673 





223 





133 





135 





82 





177 





94 





740 





442 


BP Tau 


11 


12 


130 





085 





630 





196 





311 





531 





68 





363 





84 





450 





499 


CI Tau 


7 


13 


107 





085 





874 





177 





203 





340 





72 





300 





85 


1 


635 





532 


CO Ori 


15 


10 


960 





622 


1 


473 





385 





261 





061 





56 





105 





84 





668 





425 


DF Tau 


19 


11 


983 





195 


1 


345 





292 





217 





392 





83 





386 





92 


1 


642 





629 


DC Tau 


18 


12 


207 





541 





723 





247 





341 





144 





74 





187 





91 


1 


032 





484 


DI Cep 


17 


11 


389 





059 





280 





059 





211 





190 





46 





178 





48 





973 





561 


DI Tau 


9 


12 


843 





013 





156 





064 





409 





147 





09 





440 





57 





998 





428 


DK Tau 


11 


12 


583 





178 


1 


863 





590 





316 





151 





58 





223 





90 





250 





444 


DL Tau 


9 


13 


146 





054 


1 


015 





241 





238 





189 





61 





300 





85 


1 


176 





531 


DM Tau 


6 


13 


928 





093 





591 





054 





091 





596 





79 





472 





83 


1 


285 





571 


DN Tau 


10 


12 


332 





045 





289 





081 





279 





447 





46 





223 





54 


2 


177 





542 


DR Tau 


16 


11 


660 





077 


1 


271 





371 





292 





173 





66 





229 





80 





617 





432 


EH Cep 


17 


12 


391 





248 


1 


064 





323 





303 





046 





38 





095 





63 





987 





423 


EZ Ori 


7 


11 


677 





078 





404 





211 





522 





256 





79 





236 





85 





369 





350 


GG Tau 


7 


12 


212 





039 





328 





048 





145 





696 





50 





139 





30 


1 


375 





538 


GI Tau 


8 


13 


392 





229 


1 


840 





347 





189 





231 





58 





269 





92 


1 


014 





377 


GK Tau 


8 


12 


668 





056 


1 


184 





253 





214 





180 





42 





226 





82 





747 





382 


GM Aur 


9 


12 


052 





039 





378 





103 





273 





080 





12 





369 





71 


1 


059 





500 


GW Ori 


15 


9 


878 





088 





200 





180 





901 





159 





54 





210 





79 





171 





360 


Haro 1-1 


6 


13 


355 





078 





656 





176 





268 





390 





40 





214 





52 


1 


293 





520 


Haro 1-14 


6 


14 


320 





141 





833 





170 





203 





072 





09 





156 





48 


1 


086 





433 


Haro 1-4 


5 


13 


529 





018 





291 





107 





367 





366 





19 





464 





59 





478 





609 


HI Ori 


7 


13 


347 





278 





368 





120 





326 





135 





31 





195 





70 





930 





469 


HN Tau 


6 


13 


993 





403 


1 


446 





340 





235 





055 





33 





110 





66 


1 


069 





612 


LkCa 15 


5 


12 


144 





065 





837 





289 





346 





250 





88 





188 





89 





208 





263 


LkHa 191 


5 


12 


960 





018 





178 





020 





115 





211 





21 





063 





06 





841 





421 


LkHa 209 


8 


11 


858 





026 





179 





084 





469 





241 





31 





313 





49 


1 


022 





435 


LkHa 321 


7 


12 


342 





020 





143 





029 





200 





034 





05 





118 





22 





599 





492 


LX Ori 


6 


12 


091 





039 





140 





024 





170 





031 





06 





182 





40 





800 





466 


RW AurAB 


16 


10 


516 





235 


1 


411 





289 





205 





189 





74 





229 





93 





941 





488 


RY Tau 


11 


10 


346 





334 





827 





270 





326 





033 





21 





018 





11 





956 





462 


SU Aur 


19 


9 


239 





116 





368 





230 





624 





269 





88 





240 





88 





414 





354 


SW Ori 


5 


12 


980 





190 





728 





083 





114 





181 





74 





202 





84 


1 


107 





326 


T Tau 


17 


9 


972 





062 





220 





055 





252 





240 





50 





203 





54 


1 


191 





541 


UX TauAB 


6 


11 


286 





224 


1 


487 





496 





333 





251 





96 





246 





98 





138 





342 


UY Aur 


13 


12 


560 





371 


1 


050 





365 





348 





385 





86 





314 





89 





973 





544 


V1082 Cyg 


6 


13 


142 





066 





828 





119 





144 





356 





84 





222 





78 





578 





552 


V1121 Oph 


14 


11 


430 





132 





385 





105 





272 





427 





72 





164 





62 


1 


419 





569 


V2058 Oph 


7 


12 


825 





051 





277 





083 





299 





494 





48 





342 





63 





753 





544 


V2062 Oph 


7 


12 


476 





051 





314 





218 





693 





453 





50 





341 





56 





429 





483 


V360 Ori 


5 


12 


529 





056 





152 





044 





290 





311 





37 





149 





39 





525 





496 


V521 Cyg 


17 


13 


736 





137 





732 





291 





397 





275 





61 





246 





71 





291 





289 


V625 Ori 


5 


13 


136 





112 





862 





408 





474 





385 





73 





258 





83 





364 





458 


V649 Ori 


7 


12 


078 





042 





262 





046 





176 





249 





41 





238 





62 





671 





504 


V853 Oph 


7 


13 


443 





142 





580 





203 





350 





533 





61 





544 





83 


2 


111 





596 
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